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Shock-Induced Volatile Loss from a Carbonaceous 
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ABSTRACT 



Solid recovery impact-induced volatile loss experiments on the 
Murchison C2M meteorite indicate that for an impact of a given velocity, 

H 2 0 and total volatiles are driven from the sample in the same proportion 
as present initially. Thus, the primitive surface volatile budget (volatiles 
reacted with the solid surface plus atmosphere) of a planet growing by 

accretion would have the same bulk elemental composition as the volatiles a. 

“> 

in the incident planetesimals. Incipient devolatilization of Murchison 
occurs at an initial shock pressure of about 11 GPa and complete devola- 
tilization occurs at a pressure of about 30 GPa. For the Earth, incipient 
and complete devolatilization of accreting planetesimals would occur when 
the planet reached approximately 12% and 27%, respectively, of its 
present-day radius. Thus, impact-induced devolatilization would pro- 
foundly affect the volatile distribution within the accreting planet. Prior 
to metallic core formation and internal differentiation the growing planet 
would have a very small core with the same volatile content as the 
incident material, a volatile depleted ‘mantle’, and an extremely volatile 
rich surface. In the case of the Earth, 99.4 wt% of the total incident 
volatile material would end up on or near the planetary surface. 


•now at Department of Geology, Arizona State University, Tempe, AZ 85287 
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Introduction 

Impact-induced volatile release during planetary accretion may have played an 
important role in planetary development and in atmospheric and oceanic formation [l, 2, 
3, 4]. Previous experimental results indicate that the minerals serpentine and calcite 
begin to devolatilize (i.e., lose H 2 0 and C0 2 , respectively) upon impact at initial shock 
pressures of 10 to 15 GPa and completely devolatilize at shock pressures of 30 to 40 GPa 
[5, 6, 7]. Assuming homogeneous accretion and symmetric impact, infalling 

planetesimals could have achieved impact velocities sufficient to generate shock pressures 
of these magnitudes when the Earth’s radius was about 10 to 15% and 25 to 35%, 
respectively, of its present value. Thus, extrapolation from single-mineral results sug- 
gests that impact-induced devolatilization would have been a major feature of planetary 
accretion throughout most of the growth period of an Earth-sized planet. 

We report here results of shock recovery experiments on Murchison carbonaceous 
chondrite (C2M) performed to determine the conditions for loss of H 2 0 and other vola- 
tiles upon impact from material more representative of an accreting planetesimal than 
single minerals. One reasonable form in which the volatiles of the terrestrial planets may 
have arrived is in a mineral assemblage similar to that found in carbonaceous chondrites. 
It is important to study this material directly because shock propagation, and hence 
shock-induced devolatilization, in a heterogeneous medium is complicated by local varia- 
tions in mechanical properties that can lead to localized heating, chemical reaction, and 
other types of anomalous behavior. Furthermore, the phyllosilicates that comprise the 
bulk of the matrix of Murchison include a significant proportion of Fe-rich material [8], 
whereas the serpentine shock recovery work performed to date is confined to the Mg end 
members antigorite and lizardite [5, 7, 9). Volatile-containing minerals present in 
Murchison are listed in Table 1. 



- 3 - 


Experimental 

The experimental procedure used was similar to that employed in previous solid 
recovery shock-induced devolatilization studies [5]. Murchison samples (discs 4.4 mm in 
diameter, 0.5 mm thick) were polished until the surfaces were flat and held in vented 
target assemblies that exposed the samples to the ambient atmosphere during the 
impact. The assemblies and flyer plates were made of stainless steel, type 304. The 
experiments were performed on the 20mm propellant gun at the California Institute of 
Technology. Sample shock pressure was calculated using a one-dimensional impedance 
match method [10]. The Murchison Hugoniot was assumed to be the same as that deter- 
mined by R. G. McQueen for a ‘serpentine’ of similar density (see [11]). The shock wave 
equation of state parameters used for Murchison were p 0 = 2.85 g/cm 3 , C 0 = 5.80 km/s, 
and S = 0.62 (where U B = C 0 + S u p , U s is the shock velocity, and u p is the particle 
velocity) and were obtained from a linear fit to McQueen’s U s -u p data in the range 0 
u p 3 km/s [11]. 

In the recovery assembly, the sample is sandwiched between two steel surfaces. 
During the shock event the sample is shocked first to a pressure-density point on the 
sample Hugoniot, but subsequent reverberations bring the sample to a higher pressure 
state with a density greater than that of the Hugoniot (see [12], for a description of 
reverberations in shock recovery experiments). In this work we report both the initial 
(Hugoniot) shock pressure Pj, and the peak reverberated pressure P p for each shot. In 
applying the results to models of planetary accretion, we use a formulation for devolatili- 
zation as a function of Pj (as opposed to P p ). This is because current theoretical models 
of the impact-induced devolatilization process are generally based on shock entropy cal- 
culations [3, 13, 14], and because after the initial shock, the reverberations are nearly 
isentropic. Thus, the entire entropy increase of a sample during a shock recovery experi- 
ment occurs during the first shock. Previous shock recovery experiments have often 
been reported in terms of P p [6, 15, 16]. Data used from previous studies have been 
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corrected to indicate Pj (see discussion). 

After the shock experiment, the impacted sample was recovered from the target 
assembly and analyzed for H 2 0 and total volatile content. Shock-induced H 2 0 and vola- 
tile loss were obtained by comparison with analyses of unshocked material. Total volatile 
content was measured by thermogravimetric analysis (TGA) using a Mettler Ther- 
moanalyzer 2000C. A heating rate of 5°C/min, a maximum temperature of 1200 'C, 
and a nitrogen flow rate of 20 ml/min were employed. The total volatile content of 
unshocked Murchison released over the range 200°C-1200°C was 10.09 ± 0.17 wt %. 
Water content was determined using a DuPont Moisture Evolution Analyzer with a 
maximum temperature of 1000 ° using N 2 as a carrier gas. The Moisture Analyzer was 
operated manually; the temperature was increased in increments of 50-100 'C, whenever 
the cell voltage indicated that the electrolytic cell was below the saturation level. This 
procedure was adopted to allow minimum analysis time, thereby minimizing the baseline 
correction and yielding the most precise total H 2 0 values. H 2 0 loss as a function of 
temperature under constant heating rate conditions were not obtained as a result of 
using this procedure. The H 2 0 content of unshocked Murchison used for reference in 
this work was 6.50 ± 0.20 wt % (200-1000° C). 

Results and Discussion 

Figure 1 shows representative TGA curves for shocked and unshocked Murchison. 
Unshocked Murchison meteorite loses mass over a broad temperature range, about 
260 * C-1020 * C, and the minor structure in the curve indicates that different com- 
ponents are decomposing in different temperature intervals. The TGA curves of the 
shocked material display much less structure. Mass loss begins at a higher temperature 
than for the unshocked material, 650-750 ° C, indicating that the more labile materials 
were decomposed during shock and release. No attempt was made to determine which 
components of the matrix were decomposed during the shock. 
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Total volatile loss and H s O loss results are listed in Table 2. These data are plot- 
ted in Figure 2. Incipient devolatilization occurs at an initial shock pressure Pj of about 
11 GPa, and complete devolatilization occurs at an initial shock pressure of about 30 
GPa. Aside from one apparently erroneous data point, there is no discernible difference 
between H z O loss and total volatile loss as a function of initial shock pressure. There- 
fore, there is no fractionation of H 2 0 relative to total volatiles caused by impact- induced 
devolatilization of carbonaceous chondritic material. The atmospheric (or surficial) com- 
plement of volatiles in an atmosphere created by impact induced degassing of 
planetesimals will have the same bulk elemental composition as the volatiles contained in 
the incident material, and will also be the same as that of the buried volatile component. 
Chemical reaction, both in the gas phase and with surficial materials, will change the 
actual speciation. 

It must be stressed that the above conclusion applies strictly only to those volatiles 
released in an analysis carried out in a N 2 carrier gas, at temperatures up to 1200 °C. 
Analysis under oxidizing or reducing conditions would liberate a different complement of 
volatiles. For example, graphite is stable (non-volatile) at 1200 °C in N 2 , but would be 
oxidized to C0 2 (».e., considered a volatile) in an oxidizing atmosphere. Furthermore, 
there are strong indications that the degree of impact-induced devolatilization depends 
on the ambient atmospheric conditions during both the shock experiments [7, 14, 17] and 
planetary accretion. In the experiments described here, the sample was exposed to air 
during the shock. Details of these atmospheric dependences are subjects for further 
investigations. 

Figure 2 also shows impact-induced H 2 0 loss for crystal-density serpentine [7] and 
for 20% porous serpentine [9], which approximate the range of densities of many carbo- 
naceous chondrites [18]. These results bracket the data for Murchison, consistent with 
the observation that septechlorites are the major H z O- bearing phase in Murchison. This 
also indicates that the results of devolatilization experiments on Mg end-member 
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serpentine may be extrapolated to carbonaceous chondrites. 

Laboratory shock pressures can be related to planetesimal impact velocities by 
assuming uniform, homogeneous accretion. For symmetric impact, the impact velocity 
required to generate a given shock pressure is twice the particle velocity of the shocked 
state. Minimum impact velocities onto a planetary surface are equal to the escape velo- 
city of the planet. Thus, assuming constant density, the radius of the planet is a linear 
function of impact velocity 

r = V p (8tt Gp/3)' 1 / 2 = 2 u p (8tt Gp/3)- 1 / 2 (1) 

where r is the planetary radius, V p is impact velocity, G is the gravitational constant, "p 
is the average density of the planet, and u p is the particle velocity. Figure 3 is a plot of 
shock-induced volatile loss versus fractional radius for Earth, Venus, and Mars using the 
parameterized devolatilization versus initial shock pressure results for Murchison 
presented here (solid line in Figure 2). For the Earth and Venus, impact-induced devola- 
tilization of incident planetesimals would have begun by the time the planetary radii 
were about 12% of their present-day values, and complete devolatilization would have 
occurred when the planets reached about 27% of their present day radii. For Mars, inci- 
pient and complete devolatilization of incident planetesimals would occur at about 26 
and 58%, respectively, of its present day radius. 

These results indicate that the initial volatile distribution will be extremely nonuni- 
form in an accreting planet. In the absence of thermal effects, the net result of accretion 
as outlined above is a core with the volatile content and elemental composition of the 
incident material; a volatile depleted mantle will probably form in association with a 
metallic core, the core will sink into the center of the planet, and the volatiles depleted 
from the mantle will contribute to a volatile-rich veneer at the planetary surface. The 
distribution of volatiles between the interior and surface regions depends on planetary 
radius and mean density. Integrating with respect to mass along the curves in Figure 3 
indicates that for the Earth, Venus, Mars, and a Moon-like planet, 99.4, 99.2, 93.2, and 
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42.2 mass %, respectively, of the incident volatiles end up at the planetary surface at the 
end of accretion. The possible existence of a planetary magma ocean in the later stages 
of accretion, whether caused by a runaway greenhouse effect due to the accumulating 
protoatmosphere [4] or by retention of heat deposited during impact [20], does not alter 
the conclusion that a volatile-containing core surrounded by volatile-depleted region 
exists early in planetary accretion. The existence of heat generated by core formation 
and radioactive decay, of course, profoundly affects planetary differentiation. 

Conclusions 

Impact-induced H 2 0 loss and total volatile loss from a carbonaceous chondrite have 
the same dependence on shock pressure. Thus, there is no enrichment or depletion of 
H 2 0 relative to other volatile species in gases released from a carbonaceous chondrite 
upon impact. The actual impact pressure dependence is bracketed by results for crystal 
density serpentine and 20% porous serpentine. For the Earth, incipient devolatilization 
of accreting planetesimals occurred when the planet reached 12% of its present day 
radius, and complete devolatilization occurred at 27% of its present day radius. There- 
fore, impact-induced devolatilization of accreting planetesimals played a major role in 
altering the distribution of volatiles during accretion. Prior to metallic core formation 
and internal differentiation, 99.4 mass % of the incident volatiles would be found in the 
near surface environment (including a possible surficial magma ocean), the remainder 
residing in the core region, with a volatile depleted ‘mantle* between them. 
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TABLE 1 


Volatile-containing minerals in Murchison jl9] 


Water 


Sulfur 


Serpentine* - (Mg,Fe) 3 Si 2 05(0H) 4 Hydrocarbons* 

(Septechlorites) (H/0 — l) 


Montmorillonite - 

(Na,l/2Ca).33(Al ; Mg) 2 SI 4 O 10 (OH) 2 nH 2 O 
Brucite - Mg(OH) 2 
Limonite - Fe 2 0 3 *nH 2 0 


Graphite - C 
Calcite - CaC0 3 
Whewellite - CaC 2 0 4 *H 2 0 


Sepiolite - Mg 2 Si30 8 2H 2 0 
Mg-Sulfate - MgS0 4 nH 2 0 
Gypsum - CaS0 4 2H 2 0 
Bloedite - Na 2 Mg(S0 4 )o-2H 2 0 
Whewellite - CaC 2 0 4 'H 2 0 


Disseminated S* 

Troilite - FeS 

Pentlandite - 
(Fe,Ni) 9 S 8 

**FeSO phase 

Gypsum - CaS0 4 *2H 2 0 


♦Major volatile bearing phase. 

♦♦Poorly characterized, unidentified Fe, S, O containing phase. 



Table 2. Experimental Results 


Shot 

Sample 

Projectile 

Initial 

Peak 

Total 

H 2 0 Loss, 

number 

mass, 

Velocity, 

Shock 

Shock 

Volatile 



mg 

km/s 

Pressure, 

Pressure, 

Loss, 





GPa 

GPa 

20 0 -1200 * C, 

20 0 C-1200 0 C, 






%of 

% of 






initial 

initial 

799 

19.6 

1.08 

13.3 

22.9 


14.4 



±0.03 

±0.4 

±0.7 


±3.5 

803 

19.0 

1.21 

15.1 

26.2 


40.0 



±0.02 

±0.3 

±0.5 


±5.0 

785 

19.3 

1.40 

17.3 

31.1 

40.0 

_ 



±0.10 

±14 

±2.7 

±5.0 


835 

17.6 

1.67 

21.5 

38.4 

87.4 

89.5 



±0.01 

±0.1 

±0.3 

±2.6 

±6.7 

825 

39.5 

1.84 

24.0 

43.2 

81.3 




±0.01 

±0.2 

±0.5 

±3.8 


781 

19.4 

2.06 

27.3 

49.7 

48.0 

89.0 



±0.10 

±1.5 

±3.1 

±10.0 

±5.0 



Figure Captions 


Figure 1. Representative thermogravimetric analysis (TGA) curves for unshocked and 
shocked Murchison. Mass loss increases downward. For the shocked samples, the 
impact-induced mass loss is represented by the difference between the final (1200 ° C) 
value and 100%. 100% mass loss refers to the mass loss between 200 and 1200 °C for 
unshocked Murchison. 

Figure 2. Shock-induced mass loss as a function of initial shock pressure for Murchison. 
Squares represent H 2 O loss; circles, total volatile loss. Solid line is a polynomial fit to 
the data. Dashed lines are fits to data on impact-induced dehydration of 20% porous 
serpentine [9] and crystal-density serpentine [7], 

Figure 3. Impact-induced devolatilization of Murchison as a function of fractional plane- 
tary radius and mass for the Earth, Venus, Mars, respectively. 
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